Improved Theoretical Predictions of Microlensing Rates for the 
Detection of Primordial Black Hole Dark Matter 
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ABSTRACT 

Primordial Black Holes (PBHs) remain a Dark Matter (DM) candidate of the Standard Model 
of Particle Physics. Previously, we proposed a new method of constraining the remaining PBH 
DM mass range using microlensing of stars monitored by NASA's Kepler mission. We improve 
this analysis using a more accurate treatment of the population of the Kepler source stars, their 
variability and limb-darkening. We extend the theoretically detectable PBH DM mass range 
down to 2 x 10~ 10 A/q, two orders of magnitude below current limits and one third order of 
magnitude below our previous estimate. We address how to extract the DM properties such as 
mass and spatial distribution if PBH microlensing events were detected. We correct an error in 
a well-known finite-source limb-darkening microlensing formula and also examine the effects of 
varying the light curve cadence on PBH DM detectability. We also introduce an approximation 
for estimating the predicted rate of detection per star as a function of the star's properties, 
thus allowing for selection of source stars in future missions, and extend our analysis to planned 
surveys, such as WFIRST. 

Subject headings: Black hole physics - Gravitational lensing: micro - dark matter 



1. Introduction 

Dark Matter (DM) has been the topic of exten- 
sive research and remains one of the long stand- 
ing mysteries in cosmology. Recent technologi- 
cal developments ha ve aided th e search for DM 
particle candidates ( Feng 20101) with little suc- 
cess, so increased attention is now focused on clos- 
ing the window of one of the few DM candidates 
left in the Standard Model of p article physics 
prim o rdial black holes (PBHs) (jFrampton et al 



lapse of cosmic loops, or d omain walls ([Khlopov 
2010). First proposed by Zel'dovich fc Novikov 
(|l966l ) and iHawkind (|l971f l. PBHs would form 



during the radiation-dominated era, and therefore 
would be non-baryonic, satisfying the big bang 
nucelosynthcsis limits on baryons, and would thus 
be classified as cold dark matter in agreement with 
the current paradigm. 



20101 : ICarr et alJlioioh . As opposed to regular 



black holes, PBHs are much smaller and would 
only be able to form during the early universe, 
when perturbations could collapse to form stable 
PBHs whose mass would be on the order of the 
mass of the horizon at the time of their collapse. 
PBHs can form due to density fluctuations dur- 
ing differ ent inflation scenarios , such as double 
inflation ( Frampton et al. 201dl ). as well as due 



to phase transitions in the early universe caus- 
ing a soft equation of state, bubble collisions, col- 



The disco very of Hawking evaporation (jHawkine 
197i Il975l) . led to a theoretical lower limit on 
their mass scale of 5 x 10~ 19 M Q , with any black 
holes smaller than this having evaporated by the 
current age of the Universe. There are no the- 
oretical limits on the rest of the mass range, 
however higher masses have b een progressively 
ruled out by various experiments ( Carr et "alll2010t 
lAlcock et al. 1998 ). leaving one major uncon- 
strained window left, covering 5 x 1O~ 14 M to 
2 x 10~ 8 M Q . The lower mass end of this window 
is set due to fe mtolensing observation s of gamma- 
ray bursts (see iBarnacka et al.l (|2012j ) for a recent 
analysis of the Fermi Gamma-ray Burst Monitor 
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data), while the higher mass end is constrained 
by the combi ned MACHO/ERO S r esults due to 
microlensing ( Alcock et al.l 19981 ) . In iGriest et al.l 
(|201lL hereafter Paper I), we proposed to extend 
these microlensing constraints using the NASA 
Kepler satellite, which has the capability to close 
a significant part of the window. 

The Kepler satellite is a 1 m aperture tele- 
scope with a 115 deg 2 field-o f-view in an Eart h 
trailing helioc e ntric orbit (see Koch et al. ( 2010f ): 
Bor ucki et al. I (l2010h for a description of the Ke- 
pler mission). It takes photometric measurements 
of around 150,000 stars every 30 minutes in the 
Cygnus-Lyra region. The telescope was launched 
in March 2009 and the mission has recently been 
extended to the year 2016. With planet-finding 
as its main science interest, it measures changes 
in stellar flux down to one part in a thousand or 
lower. This proves very beneficial in microlensing 
searches as well, where flux magnification is mea- 
sured in the stellar light curves. 

Here we fill in details and expand on our previ- 
ous analysis of this exciting possibility, analyzing 
the importance of limb-darkening on tightening 
our constraints, as well as calculating the proba- 
bility for various lens parameters if any events are 
detected. We also improve our numerical estimate 
over the previous use of a 5000 star sample, by 
including all the third quarter Kepler stars being 
monitored (^156,000 stars), and project this as 
representative of the majority of the mission. We 
introduce new notation for the limb-darkened mi- 
crolensing curves, which will be used to fit any fu- 
ture observed events. We correct an error in a well- 
known finite-source limb-darkening microlensing 
formula and derive a limb-darkened microlensing 
detection efficiency in our framework. Finally, we 
show that the PBH DM mass window can be ex- 
tended further to lower masses using this improved 
analysis, and provide an approximation which can 
be used to predict microlensing rates in future sur- 
veys, such as the planned WFIRST space mission. 

2. Formulas 

2.1. Point-Source Microlensing 



Paczvnskil ([19860 first proposed microlensing as 
a way to search for DM in the Milky Way. In doing 
so, he assumed a point-source, point-lens model, 
in which a lens, such as a PBH, would cause a 



magnification of this source when crossing in front 
of it, described by 

„,2 , 9 



-4„ 



(1) 



where u = b/rs and b is the impact parameter of 
the lens, that is the transverse distance between 
the lens and the line-of-sight to the source. The 
Einstein ring radius te is given by 



?"E 



4GMLx(l ~ x) 
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where x is the ratio of the lens distance to the 
source distance, L is the distance to the source 
star, and M is the mass of the lens. As the PBH 
passes in front of the star, the amplitude becomes 
time-dependent, A{t) = A[u(t)}, and 

1/2 



u(t) 



2(t - to) 



-I 2 



(3) 



Here to is the time of the peak magnification, 
u m in = u(to), and £e = 2rE/ft is the time for 
a lens to cross the Einstein ring with a velocity vt 
transverse to the line of sight. 

This is the standard point-source limit, in which 
the lens produces a 34 percent magnification when 
it is within one Einstein radius of the projected 
source star. This approximation is valid for a 
source that is much smaller than the Einstein ra- 
dius and is not directly aligned with the lens. How- 
ever, for the relatively nearby Kepler source stars 
and the relative low mass PBHs, the projected star 
radius needs to be taken into account. 

2.2. Finite-Source Microlensing 

When the projection of the radius of the star 
is comparable to the lens impact parameter, one 
needs to take into account finite-source effects on 
the detectability of events. The projected star ra- 
dius is given by 
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(4) 



ve(x) ' 

where i?„ is the radius of the source star. For a 
consta nt surface brightne ss, equation [1] now be- 
comes (|Witt fc Maolll994 eqns 9-11) 
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for u = {/*, and 
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F, E, and II are elliptic integrals of the first, sec- 
ond, and third kind. As opposed to the point- 
source approximation, there is now a maximum 
amplitude for the m agnification, which is equal to 
(jWitt fc Maolll994l eqn 13) 



(9) 



The duration of the event i e vent will now be 
the time during which the event is detectable, 
starting when the lens crosses the threshold im- 
pact parameter Wthresh at which the microlensing 
light curve magnification reaches the minimum de- 
tectable threshold Athresh = A(u t hrcsh, U*). The 
duration of the event is then described by 
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As the ratio of the distance to the lens to the 
distance to the star, x, approaches 1, the pro- 
jected star radius in terms of the Einstein radius, 
U*, approaches oo, thereby suppressing the max- 
imum magnification A maXft . Therefore there is 
some x max beyond which A maxts is lower than the 
detectable magnification, A t hresh and no events are 
detected. This effect decreases the detection effi- 
ciency of the PBHs. However, the duration of the 
event increases, since the PBH does not have to 
be one Einstein radius away from the center, but 
from the edge of the projected star, for the begin- 
ning of a microlensing event. This effect increases 



the detection efficiency As seen in Paper I, this 
dominates the detectability of events, increasing 
the number of expected microlensing events in the 
Kepler light curves. This finite-source model, how- 
ever, assumes a constant brightness of the star, 
and does not take into account limb-darkening of 
the source star. In this paper we extend our anal- 
ysis of the detectability of events to include this 
more physical model of the stars being lensed. 

2.3. Finite-Source Microlensing with Limb- 
Darkening 

For microlensing of nearby Kepler stars, where 
the Einstein radii of detectable PBHs is very small, 
limb-darkening is anticipated to be an impor- 
tant effect on the rate of detection. The effect 
is such that the star appears to be brighter to- 
wards the center, producing a more concentrated 
source brightness, mimicking a model in between 
the point-source and the finite-source approxima- 
tions. The limb-darkening profile we use to study 
this is the linear limb-darkening described by Witt 
and Mao, projecte d into the lens plane where 
dWitt fc Madll994L eqn 15) 



h(Ul) = l-u x + u x yjl- (C/;/C/*) 2 . (11) 

Here, t/„ is now the distance from the center of 
the projected star in terms of the Einsten radius 
and u\ is the linear limb-darkening coefficient. 
Witt and Mao calculated a limb-darkened profile 
numerically using a weighted surface brightness 
(|Witt fc Madll994l eqn 16), 



A nmh (u, [/*) 



•u. 



2TrI h (U*)dU t 



dUL 



(12) 



The integrand in the second integral of equa- 
tion [T^] has a peak at u = U m , which causes some 
problems with convergence when integrating nu- 
merically. Wit t and Mao provide d another form 
of equation [TJ (jWitt fc Maolll994 first half of eqn 
16) but recommended against its use due to, they 
said, the presence of a singularity. We find in fact, 
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that there is an error in their equation, due to the 
treatment of the integral boundaries. Integrating 
equation [12] by parts we find the correct magnifi- 
cation for a linear limb-darkened profile to be 



Aimb{u,u*) = (i - ux/zy 



A fs (u,U,)TrU!(l~u x ) 



u. 



+ (7ra A /C/*) / Ms{u, ~ z 2 Wt - z 2 )dz 
Jo 
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where 

A fs (y)iry 2 = 2y + (1 + y 2 ) ( - + ; nvsm 
for u = y, and 



u;-i 

■2 "U 2 + l / 
(14) 



2 2{u-y) 2 l + y 2 
A&(u,y)iry = — — ; — { = U[-,n,k 



(15) 



(u + y) v /4+(u-y) 2 " V2 
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for u ^ y. Here y = f7* or y = ^/U 2 — z 2 , where 
appropriate in equation 1131 Use of this equation 
removes the problem of numerical convergence in 
equation [12] 

As stated above, the limb- darkening of the 
source star produces a more concentrated source 
brightness, thereby changing the shape of the mi- 
crolensing light curve. There will therefore be 
a higher maximum magnification than that pro- 
duced by a pure finite-source light curve. Agoll 
(|2002h calculated this maximum amplification 
for a quadratic limb-darkening model. Here we 
present the result for the linear limb-darkening 
profile using equation 1131 when u = 0: 



A„ 



1 - ua/3 
X i 1 - ux + (2u x /3U 2 ) 



(2 + U 2 )E [U 2 /(A + U 2 ) 
2K [U 2 /(4 + U 2 )] 
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where K and E are the complete elliptic integrals 
of the first and second kind. 

Therefore, there is a new, higher x max below 
which the magnification is detected, allowing more 
PBHs to be observed. At the same time, since the 
brightness of the star is more concentrated, the 
impact parameter for which an event is detected 
will be closer to the projected center of the star, 
decreasing the duration of the event. Since this is 
an effect in between the point-source and finite- 
source model, there should be an overall decrease 
in the number of expected events in the Kepler 
data. 

3. Effect of Limb Darkening on the Nu- 
merical Estimate of Expected Number 
of Events 

In Paper I we calculated the expected number 
of PBH microlensing events in the Kepler data us- 
ing a subset of 5000 stars of the third quarter light 
curves. Here we make a more accurate estimate of 
this number. We first extend the calculation to in- 
clude the full set of the third quarter light curves 
being monitored, reflecting a more accurate sam- 
ple of the stars. In addition, we scale the obser- 
vation time to 7.5 years to include the extended 
Kepler mission, as well as assume a more accu- 
rate number of variable stars in the Kepler mis- 
sion, with 25% of the observed dwarf stars and 
95% of t he observed giant stars assumed to be 
variable (|Ciardi et al.ll201l[ ). This results in a to- 
tal of 780,000 star-years being observed, twice the 
amount assumed in Paper I, however now the dis- 
tribution is assumed as consisting of mainly dwarf 
stars with many of the giant stars being unus- 
able. The previous estimate of the number of 
expected PBH microlensing events cannot there- 
fore be naively multiplied by two. As in Paper 
I we require 4 sequential measurements 3-sigma 
above average, equivalent to a microlensing event 
of a minimum 2 hours duration. Calculating the 
stellar distances and the magnification thresholds 
^thresh (below which an event would not be de- 
tected) for each star as described in Paper I, we 
recalculate the number of expected events using 
equation 1 of Paper I, for a finite-source model of 
microlensing, with this new distribution of stars. 
This is plotted in Figure [T] as the solid line. The 
shape and peak of the number of events for each 
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Fig. 1. — Top panel: Expected number of events, 
iVe X p, scaled to 780,000 star-years of Kepler ob- 
serving time for a finite-source microlensing model 
with no limb-darkening (solid line) and with limb- 
darkening (dashed line). We assume 4 sequential 
measurements with 3-sigma above average mea- 
surements of the flux. The dash-dotted line repre- 
sents N cxp from our previous estimate of Paper I. 
The horizontal line shows the 95% confidence level 
limit if no events are detected. Bottom panel: The 
potential 95% confidence level exclusion of PBH 
dark matter. The area above the solid line for 
the finite-source model (dashed line when limb- 
darkening is included) would be ruled out if no 
events are detected in the 7.5 year Kepler obser- 
vation lifetime. The dash-dotted line represents 
the previous estimate of Paper I, while the dotted 
line represents the current limits from the com- 
bined MACHO /ER OS LMC microlensing surveys 



( Alcock et al. 19981 ). The horizontal line depicts a 



DM halo consisting entirely of PBHs. 



PBH mass has changed. Discarding 95% of the 
giant stars as variable results in the curve being 
dominated by the dwarf stars' sensitivity. While 
the giant stars' peak PBH DM mass sensitivity is 
at around 10 _8 M Q , which dominated the earlier 
curve in Paper I, the dwarf stars' sensitivity is bet- 
ter at the lower PBH masses, peaking at around 
5 x 1O _9 M0, since a higher fraction of the stellar 
intensity is magnified. The larger amount of star- 
years, dominated by monitoring of dwarf stars, in- 
creases the sensitivity to lower mass PBHs. This 



more accurate estimate shows a potential closing 
of the PBH DM window down to 2 x 10- 10 M Q , 
compared to the 5 x 10 -10 previously estimated. 

In addition to these calculations, we now con- 
sider the effect of limb-darkening on the predicted 
number of detectable events. As seen in the pre- 
vious section, limb-darkening increases A max , but 
decreases the duration of the event, since the to- 
tal stellar flux remains the same. A max determines 
£max) the distance to which a PBH would be de- 
tectable, with 1 being the maximum value. A 
higher A max would naively increase this value, in- 
creasing the range of masses that PBHs would be 
detectable, however, calculating x niax for a typ- 
ical Kepler star (with radii between 0.9i?o and 
1.5i?0), without limb-darkening, one can see that 
i mEU[ is already approaching the maximum value 
of 1. Therefore, increasing this further when limb- 
darkening is added, does not produce much of an 
effect on the total number of events that could be 
detected. The only other effect then is to reduce 
the possible duration of an event, therefore de- 
creasing the total expected number. In Paper I, we 
showed that the naive point source optical depth 
(the total number of PBHs inside a microlens- 
ing tube as defined in iGriestl (Il99lh l is propor- 
tional to ^thresh (^ ne detectable impact parameter 
value) and thus replacing this with the projected 
star radius, changes the optical depth by a factor 
of U*. Extending this naive calculation to limb- 
darkening, we can see that any effect on the pro- 
jected star radius will directly impact the optical 
depth in quadrature. Therefore, limb-darkening, 
which effectively reduces the radius of the star, 
could drastically limit how far we can extend the 
current PBH mass range, since lower PBH masses 
would not be as detectable. 

Here we explore this effect by repeating the 
above calculation for the expected number of PBH 
microlensing events, however now including the 
limb-darkening effect for each star. In order to 
do this, we used the linear limb-darkening model 
as described in Section [2~3l calculatin g the linear 
limb-darkening coefficients using the (|Sindl20inh 
model grid to find Kepler limb-darkening coef- 
ficients as a function of the effective tempera- 
ture, surface gravity, and metallicity of each of 
the 150,000 Kepler source stars. This enables 
us to calculate a new detectable magnification 
threshold A t hresh using the magnification formula 
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in equation [13] for linear limb-darkening. The 
number of expected events including this linear 
limb-darkening effect is plotted in Figure [1] as the 
dashed line. 

Surprisingly, Figure Q] does not show a signif- 
icant effect of limb-darkening on the number of 
expected events. This can be understood as being 
due to the extreme precision of Kepler. Since Ke- 
pler light curves allow a magnification threshold 
of thresh = 1-001 or lower to be set, the median 
Kepler star (with a small stellar radius), will al- 
low for a detection as soon as a PBH is within 
one Einstein radius of the projected star radius, 
where the limb-darkening does not play much of 
an effect yet on the light curve. On stars with 
bigger stellar radii, as well as models with smaller 
PBH masses, this limb-darkening will have some 
effect, since the magnification reaches ^thresh only 
when the PBH's impact parameter is inside the 
projected star radius, where the limb-darkening 
plays a role. However, even then, the duration is 
shortened by just a small amount. In Figure [2] we 
plot the percent change between the finite-source 
model with and without limb-darkening and we 
in fact see that there is a decrease of up to 17% 
for the lower lens masses. Moving towards higher 
PBH masses, almost all the Kepler stars which 
have radii on the order of a solar radius, will be 
able to have a detectable event within one Ein- 
stein radius of their projected radius, therefore, 
making the limb-darkening effects negligible. We 
found that if we did not discard 95% of the gi- 
ant stars as variable, we would see an increase 
in the overall number of expected events due to 
limb-darkening at higher mass ranges. A max is 
significantly lower for stars with large radii, and 
therefore the detectable value of x will be close to 
0. We found that introducing limb-darkening in 
giant stars increases this ^4 max , increasing the pos- 
sible detectable range of x, drastically increasing 
the number of detectable events. Also, the higher 
the mass of the PBH, the lower the projected star 
radius, and the more pronounced this effect will 
be, increasing the x ratio by a higher amount when 
introducing limb-darkening. Therefore the biggest 
effect of limb-darkening on the number of events 
seems to be in these large-radii stars and on the 
lower lens masses. The effect of limb-darkening 
will also be crucial in the fitting of potential PBH 
events as well as in calculating the experimental 



detection efficiency. 
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Fig. 2. — Percent change in the number of ex- 
pected microlensing events due to limb-darkening 
for 780,000 star-years of Kepler observing time, 
defined as iVe X p with limb darkening divided by 
-ZVexp without limb-darkening included. The hori- 
zontal line depicts 100%, representing no change, 
for comparison. 



4. Detection Efficiency 

In order to turn either detection or non- 
detection of microlensing into a statement about 
Galactic dark matter, one must calculate the ef- 
ficiency at which the experiment finds PBH mi- 
crolensing. This is done by performing a Monte 
Carlo simulation with randomly generated mi- 
crolensing events and using the same selection cri- 
teria used to select microlensing light curves. One 
thus calculates the detection efficiency e, which 
is just the fraction of simulated events recovered. 
Using e one can calculate the expected number of 
detectable events as 



N cxp = 



JO 

dT 



e(x,u m i n ,v t ) 



dxdUmindVi 



■dvt du m i n dx, 



where 



dT 



dxdUmi-ndVt 



4r E (x)L 



(17) 
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where v c « 220 km/s is the halo circular velocity, 
p « 0.3 GeV cm~ 3 is the local dark matter density, 
and vt is the transverse lens velocity. Here we are 
making the approximation that the DM density is 
constant between the Earth and the source stars, 
valid for the relatively nearby Kepler stars, and 
using the fact that the Kepler field is nearly in the 
direction of the Sun's motion around the galaxy 
(see Paper I). 

Without limb-darkening, e is a function of only 
x and vt, and the integral over u m i n can be per- 
formed since that distribution is well-known due 
to a uniform stellar intensity across the projected 
star radius. This is not true when limb-darkening 
is included, since a lower u min produces a higher 
amplitude in a limb-darkened light curve as well 
as decreases the duration of an event. The effect 
of limb-darkening on the efficiency is calculated 
by adding the parameter it m i n to the Monte Carlo 
simulation. Thus the two competing effects can 
both influence this probability of detection. 

5. Detectable Parameters 

Also of interest is what we can say about a po- 
tential PBH if we do detect microlensing events. 
Assuming the measured light curve parameters are 
U* and tevent, one can then calculate the mass 
probability and the distance probability of the lens 
in terms of these measured parameters. Perform- 
ing a change of variables in terms of the observable 
parameters, the mass and distance likelihood func- 
tions can be derived. We find the mass likelihood 
function to be 



dT pv 2 c 2 R*U, ( c 2 Rl r ., . 

1 +U 2 A (3 2 g{(3), 



dt QVCnt dU* 2G M 2 \AGML 
where 



and 



±u th r Csh (U*,u x ) 2 RiU 2 ( c 2 Rl t u2 



t 2 V 2 
event c 
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(20) 



g(fi)= / dyy 3 /*(l-y)- 1 ' 2 e-to. (21) 
Jo 

In Figure [3] we plot this likelihood for £ 0VC nt = 2 
hours, R* = IRq, L = 0.73 kpc, u\ = 0.61, and 



^thresh = ^Urob (""thresh, U*) = 1.0007, values typi- 
cal of a median star in the Kepler field being mon- 
itored. The different curves represent the range of 
different values of f7* that could be measured for 
such a star undergoing a microlensing event. The 
distributions are normalized to have unit area un- 
der each curve so that each curve can be thought of 
as a probability density, that is for a measured du- 
ration of 2 hours and the given U* , the curves give 
the relative likelihood of the event being caused 
by a PBH of mass shown on the abscissa. The 
rise of each curve at low mass is dominated by the 
/3 2 term, which increases with increasing U*. The 
event duration is proportional to this projected 
star radius £/*, the PBH Einstein radius, rg, and 
inversely proportional to the PBH transverse ve- 
locity, vt ■ For a lower mass PBH only small trans- 
verse velocities will give rise to events longer than 
2 hours. If £/* is decreased, then for the lower 
mass PBHs the transverse velocities have to be 
smaller, with the number of possible events for 
these mass ranges approaching zero. This is the 
reason that the distributions tend to zero for the 
smaller values of U* for the smaller PBH masses in 
Figure [3] On the other hand, the decrease in each 
curve at larger PBH masses is caused by the lower 
number density of PBHs as their mass increases. 
The bigger the PBH lenses, the lower the number 
density needed to explain the local dark matter 
density. This in turn corresponds to less potential 
microlensing events. If microlensing events were 
to be observed in the Kepler data, we would then 
be able to, using distributions such as these, es- 
timate the mass of the PBH making up the DM. 
The product of these likelihood functions would 
give us an estimate of the PBH DM mass range. 
We could also use these distributions to exclude 
some range of masses that the microlensing lenses 
could represent by measuring their i CVO nt and U* 
parameters. 

In a similar way, we can calculate the lens dis- 
tance likelihood function, with a change of vari- 
ables in Equations [19] - [21] 



dr v 2 T 2 

= &Gp^ w U4l~x) 2 f3 2 g(/3), (22) 



where 
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to determine the distances to the lens. 
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Fig. 3. — Lens mass relative likelihood for t cvcnt = 
2 hours, = li? Q , L = 0.73 kpc, u x = 0.61, 
and ^4iimb(uthrosh, U*) = 1.0007 representing a me- 
dian star in the Kepler field. The range of values 
for U* = 10, 30, and 50 represents the range that 
could be measured in such a star. The curves here 
are normalized to have unit area under each curve, 
so that each curve represents the probability of 
each mass given the measured lightcurve parame- 
ters. 
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Fig. 4. — Lens distance relative likelihood for 
Movent = 2 hours, i?» = li?©, u\ — 0.61, and 
^4iimb(wthrosh, U*) = 1.0007 representing a median 
star in the Kepler field. The range of values for 
£/* = 10,30, and 50 represents the range that 
could be measured in such a star. The curves here 
are normalized to have unit area under each curve. 
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These distribution functions are plotted nor- 
malized to unit area under each curve in Figure 
IH We see that the distance to the lens is not very 
dependent on the value of [/* measured. The dis- 
tance probability distribution is dominated by the 
transverse velocities that are detectable at each 
distance to the lens (closer lenses have to be trav- 
eling slower in order to be detected for a given 
measured t cvcnt ). If a microlensing event was to 
be detected, then the probability distribution for 
the distance to the lens can be plotted for the par- 
ticular stellar radius being monitored, just like in 
Figure 2) A bigger i?» will yield a distribution 
centered at lower x, while a smaller R* will shift 
the distribution to higher x. We can then narrow 
down the most likely position of the lens based on 
the likelihood of detection at each distance for the 
particular star at which the event occurred. Over- 
all, though, Figure S] shows that it will not be easy 



6. Limitations 

We turn now to considerations of what could 
be improved in upcoming missions, and what the- 
oretical limitations the Kepler satellite has on de- 
tecting PBHs. As mentioned above, for the lower 
mass range of PBHs, the higher-velocity objects 
would not be detectable, because they result in 
magnifications that last too short a time. In Fig- 
ure [5] we plot the maximum detectable velocity for 
a given PBH mass for a Kepler star with median 
parameter values of distance, radius, and j4thresh, 
calculated from the third quarter stars monitored. 
Throughout this section, by velocity we mean the 
velocity of the lens relative to the Earth-source 
line-of-sight in the direction perpendicular to the 
line-of-sight. For a median type Kepler star, the 
v c value of 220 km/s is on average detectable for 
masses above 2 x 10~ 7 solar masses. The lower- 
mass range events would be dominated by the pro- 
jected radius of the star and therefore their max- 
imum velocity curve approaches a constant value, 
but at the upper-mass range, the Einstein radii 



8 



of the PBHs become important in detecting the 
event, and therefore more of the velocity distribu- 
tion is detectable. 



10000 




Fig. 5. — Maximum detectable transverse veloc- 
ity, i? maX ) of PBHs, for a median type Kepler 
star with i?* = IRq, L = 0.73 kpc, u\ = 0.61, 
and j4iimb(wthrcsh, U*) = 1-0007, for a minimum 

tevent = 2 hours. 

In order to understand the limits on this de- 
tectable velocity, we calculated the optical depth, 
which determines how many measurements are 
needed for any given PBH mass, and is defined 
as the number of PBHs inside a detectable mi- 
crolensing tube (as defined in Griest 1991, with 
the addition of an x max cutoff for the finite-source 
model). We therefore arrive at the equation t = 
tt( P /M)L f* majc u 2 thrcsh (x)rl(x)dx. In Figure © we 
plot the optical depth, averaged over all the third 
quarter Kepler source stars. The overall ampli- 
tude of the curve is set by the average distance to 
the stars being monitored, whereas the shape of 
the curve at lower-mass PBHs is mostly governed 
by x max , the detectable distance to the lens. This 
is set by the photometric accuracy of Kepler and 
cannot be changed. 

Another limitation is related to our detection of 
events that last 2 hours or more, which is set by the 
Kepler satellite cadence. In Figure [7] we plot the 
average event duration, defined as t eV ent = r /r, 
averaged over all the third quarter Kepler stars 
being monitored. We can see that the average 
event duration is about 2 hours at around 10 -8 
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Fig. 6. — Average optical depth for the third quar- 
ter Kepler source stars being monitored. The hor- 
izontal line depicts the average optical depth if a 
point-source microlensing model is used for com- 
parison. 



solar masses. With the 2 hour selection criteria, 
we are therefore not able to detect most events 
for the smaller mass PBHs. Figure [7] shows that if 
events as short as 0.1 hours could be detected, one 
might be able to detect PBHs of masses down to 
1O~ 1O M0 or below. This could be improved upon 
by decreasing the Kepler satellite cadence when 
monitoring the stars. 

Our theoretical detectable limits are due to 
both the threshold detection limit of Kepler and 
the cadence. The first cannot be changed, however 
we address improvements to the second limitation 
in the next section. 

7. Future Planned Missions 

Kepler's extended mission of four additional 
years helps increase its sensitivity to lower-mass 
DM PBHs, as seen in Section [3] Here we address 
whether any possible additional measures could be 
undertaken during this extended mission to fur- 
ther increase this sensitivity. As seen in the pre- 
vious section, in the analysis of the average mi- 
crolensing event duration, cadence plays a huge 
role in these measurements. The Kepler camera 
actually takes one image every minute. Due to 
communication bandwidth, for most stars, the Ke- 
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for 780,000 star-years with the same assumed stel- 
lar variability as in Section [3] Although requir- 
ing short cadence for all the Kepler stars for the 
full mission is not achievable, this demonstrates 
that reducing the cadence time on the current 
source stars with the current light curve precision, 
would give only a modest increase in the sensi- 
tivity to lower-mass PBHs. However, if Kepler's 
aperture was large enough to maintain the cur- 
rent signal/noise at a one minute cadence, great 
improvement would be possible. In this we note 
that certain stars are far more valuable than oth- 
ers for detecting microlensing. 



Fig. 7. — Average event duration, < t cvcnt >, for 
the third quarter Kepler source stars being moni- 
tored. The horizontal line depicts t cvcn t = 2 hours, 
which represents the minimum event duration re- 
quired for a microlensing detection in the Kepler 
data. 



pier team adds up the one minute exposures into a 
30 minute exposure before transmitting the data 
to Earth. These are called "long cadence" light 
curves. However, for a selectable subset of stars, 
the entire one minute cadence ("short cadence") 
light curves are transmitted. Thus we wish to in- 
vestigate the value for microlensing of the Kepler 
team returning short cadence data on a subset of 
stars. If we naively decrease the monitoring ca- 
dence to 1 min for all the Kepler source stars, 
for the full 780,000 star-years, the lower-mass DM 
PBH sensitivity increases by an order of magni- 
tude, down to 2 x 10~ u Af Q . Thus, while not 
possible due to bandwidth limitations, naively re- 
turning short cadence light curves for all the stars, 
would allow the exploration of one additional fac- 
tor of 10 in the allowed PBH mass range. This 
could be an exciting possibility. However, when 
the cadence is decreased, the Poisson average error 
in each flux measurement increases. Thus with the 
set light gathering power of the Kepler telescope, a 
shorter cadence is offset by a larger Athreuh< in or_ 
der to investigate this trade-off we redid our anal- 
ysis assuming a 1 min cadence, but reducing the 
signal/noise for each flux measurement appropri- 
ately. Figure [8] shows the results of this analysis 



1 min cadence 
30 min cadence 




log 10 (M/M o ) 



Fig. 8. — Total number of expected events for 
780,000 star-years for 1 min and 30 min Kepler ca- 
dence monitoring times. The horizontal line shows 
the 95% confidence level limit if no events are de- 
tected. 



We therefore address the characteristics of 
source stars and cadence times that would op- 
timize this sensitivity in future missions similar to 
Kepler. Since 95% of the giant stars are seen to be 
variable, we focus on the dwarf stars, which would 
yield a higher number of less variable light curves 
for the same number of stars monitored. To inves- 
tigate an optimal selection of stars for monitoring 
we calculated the expected number of microlens- 
ing events for each non-variable dwarf star in the 
third quarter Kepler data. We plotted this per- 
star-A^oxp versus other stellar parameters to see 
which correlated well with higher -ZV cxp . In Figure 
[51 we show the best such correlation, N cxp vs. T e g, 
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where stars with higher T e g are much more likely 
to return a detectable microlensing event. This 
figure shows that source stars with T e g > 8000A" 
are roughly 100,000 times more valuable for mon- 
itoring than stars with T e g < 4000 A. Thus, short 
cadence monitoring of a handful of carefully se- 
lected stars should be as valuable as large numbers 
of typical stars. The bandwidth problem, there- 
fore, might be solved by long cadence monitoring 
of most stars but short cadence monitoring of a 
small sample. 

Why does T c g correlate so well with the ex- 
pected number of lensing events? We see that 
for a point-source microlensing model detectabil- 
ity correlates well with luminosity; the more lu- 
minous stars give higher rates of predicted detec- 
tions when the projected star radii and the lens 
Einstein radii are comparable. However, for these 
small PBH masses, the Einstein radii are much 
smaller than the projected star radii and the lu- 
minosity per unit area is more important, giving 
the upward trend with effective temperature. 
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Fig. 9. — Expected number of events per star-year 
for a PBH mass of 10~ 10 Mq with a 1 min cadence, 
for each non- variable dwarf star being monitored 
in the third quarter of the Kepler data, plotted 
with respect to T c g. 



In order to understand better the factors that 
influence the number of predicted events for a 
given star, we can make the following approx- 
imations. For light curves with 1% or bet- 
ter precision, as seen in Kepler, we can ap- 



proximate Mthrcsh ~ U*- Also, for R^/Rq < 
0.57 (teveat/lhr) (v c /220km/s) /x max , we can ap- 
proximate (3 < 1 and /3 2 g(/3) « (3/8)n(3 2 . Using 
these approximations in equation 1 of Paper I, 
with £ max < 1 for M = 10 _ 10 Mq, we arrive at 
the predicted rate of detection for a given star 



T « 409.6tt 



G 5 M 4 p [/* max L 6 



(24) 



v 2 t 3 ■ R e ' 
where C/» max corresponds to the maximum [7* de- 
tectable for a given A t hresh of a star. Calculating 
this for the appropriate p and v c , we arrive at 
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(25) 



Using the third quarter Kepler stars, we plot this 
approximation as a straight line in Figure [TU] along 
with the actual rates calculated using our integral 
formulas. The approximation works very well for 
these low-mass PBHs and demonstrates how the 
stellar characteristics come into play in this cal- 
culation. As seen in Equation [55l a high stellar 
distance-to-radius ratio is important. The depen- 
dence of the rate on the effective temperature is 
also readily explained by the fact that the dis- 
tance is directly calculated from this value us- 
ing L = 1.19 x lO- 3 i?»(T cff /T ) 2 lO°- 2 ( y - Av+B C ) 
where V is the apparent visual magnitude, Ay 
the extinction parameter, and B.C. the bolomet- 
ric correction as in Paper I. We can see that the 
rate will be related to the effective temperature 
as T a TgJ. Also, decreasing the cadence will 
have a cubic effect on the rate expected, while 
maintaining a low A t hrcsh is also important, as 
expected. We can therefore predict, for exam- 
ple, that for a Kepler type mission, with 30 min 
cadence, in order to push the PBH mass limits 
down to 5 x 10 _11 Mq, one would have to observe 
160,000 dwarf stars with a (Z/lkpc) frac- 
tion of 3.15 or higher, while maintaining the 0.1% 
lightcurve precision. It seems therefore that push- 
ing to smaller PBH masses in order to close the 
remaining PBH DM mass window will be diffi- 
cult. However, if events were detected, a survey 
pointing towards or away from the Galactic cen- 
ter would provide us with more information about 
the DM distribution. We therefore address the 
WFIRST mission next. 
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Fig. 10. — Expected number of events per star- 
year for a PBH mass of 10~ 10 M Q with a 1 min 
cadence, for each non-variable dwarf star being 
monitored in the third quarter of the Kepler data, 
plotted with respect to the stellar variables gov- 
erning this number. The straight line represents 
the approximation in Equation 25. 



7.1. WFIRST 

NASA's Wide-Field Infrared Survey Telescope 
(WFIRST), is a proposed next-generation space 
observatory being designed to search for Dark En- 
ergy and extrasolar planets ( Green et al.ll2011 ). It 
is designed to perform Dark Energy measurements 
using Baryon Acoustic Oscillations, Type la Su- 
pernovae, and Weak Lensing. In addition, it will 
complement the Kepler mission with its microlens- 
ing search for extrasolar planets, targeting stars 
towards the galactic bulge. Here we address its 
value for PBH DM limits or characterization. The 
preliminary specifications are to monitor 2 x 10 8 
stars with a cadence of 15 min and a 1% pho- 

2010I ). It is es- 



tometry precision (jBennett et al 



pecially exciting, as it will monitor the center of 
the Galaxy, and therefore could potentially pro- 
vide insights into the DM distribution of the Milky 
Way. It will add to the existi ng DM dynamica l 
constraints due to microlensing ( Iocco et al1l201ll) . 
Here we provide a preliminary calculation for the 
number of expected events, if the stars being mon- 
itored are similar to those of Kepler, as well as 
assuming a simple DM distribution of the form 



with a — 5 kpc. Following our analysis for the 
Kepler mission, we find the rate of detection, 



r/r 



a 2 + rl 2 



dx- 



+ L 2 (l - x) 



(27) 



with all parameters as defined in Paper I. In Fig- 
ure [TT] we plot the number of expected events for 
one year of monitoring of Kepler-type stars (with 
ss li?©) toward the galactic bulge (L — 8 kpc), 
with 25% of stars assumed to be variable, and ig- 
noring transverse velocity. Alongside, we plot the 
approximation given in Equation [51] as a dashed 
line, demonstrating its usefulness for predicting 
microlensing rates for PBH lenses with masses less 
than 10 -10 Afo, as stated in the previous section. 




■11 -10 -9 -8 -7 -6 
log 10 (M/M s ) 



-5 -4 



Fig. 11. — Total number of expected events to- 
wards the galactic bulge for a mission such as 
WFIRST. The dashed line represents the approx- 
imation of Equation 25, appropriate for PBH lens 
masses below 1O -1O M0. The horizontal line shows 
the 95% confidence level limit if no events are de- 
tected. The dotted line represents the predicted 
Kepler microlensing limits as shown in Figure 1, 
for comparison. 



This is a preliminary analysis, however, it 
demonstrates that WFIRST will complement Ke- 
pler in the same PBH mass range, providing a 
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greater number of events, and exploring an addi- 
tional order of magnitude of the PBH DM win- 
dow. If the stars being monitored are larger than 
the Kepler type stars, the number of expected 
events will decrease, as predicted by Equation 1251 
Also, the transverse velocity of the source stars 
will have a sizable effect, as well as lensing due 
to other stars. If a PBHs are detected by Kepler, 
WFIRST will be able to study the DM character- 
istics, such as velocity and spatial distributions. 

8. Conclusions 

In this theoretical paper, we re-addressed the 
possibility of detecting or ruling out PBH DM us- 
ing the microlensing of Kepler source stars. We 
introduced a more proper treatment of the popu- 
lation of source stars and their variability, includ- 
ing a finite-source microlensing framework which 
includes limb-darkening. Using this analysis, we 
found that the PBH DM mass constraints could 
be extended down to 2 x 10~ 10 Af Q , two orders 
of magnitude below current limits and a third of 
a magnitude lower than our previous work. As 
mentioned, a proper Monte Carlo simulation will 
be needed to fully understand the experimental 
detection efficiency. We provide a limb-darkened 
microlensing framework which will be used to fit 
potential Kepler light curves, and which will help 
to distinguish between PBH's and stellar flares, 
the main source of experimental systematic error. 
Our analysis leaves us optimistic in being able to 
provide a probability for the masses of the lenses 
if any microlensing events are found, and therefore 
characterizing the DM and its epoch of formation. 
We found a very strong correlation between the 
rate of PBH detection for a given star and its 
effective temperature, providing an approximate 
expression for this rate for these low-mass PBHs. 
This should help in selecting stars to be monitored 
in future microlensing experiments. Using this ap- 
proximation, it can be seen that fully closing the 
remaining PBH DM window using a microlensing 
method will be difficult, however sensitivity could 
potentially be improved by another order of mag- 
nitude using future planned missions, such as the 
WFIRST survey towards the galactic bulge. More 
analysis for this case is needed. 
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